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Abstract 

The scaling rules of the invariant yields and cross sections for hard scattering 
processes in proton-nucleus ipA) and nucleus-nucleus {AB) reactions at LHC 
energies relative to those of nucleon-nucleon A^A^ (isospin averaged pp) col- 
^ ■ lisions are reviewed within the Glauber geometrical formalism. The number 

CN . of binary A A inelastic collisions for different centrality classes in p-(-Pb and 

O \ Pb-i-Pb colhsions at ^snn - 8-8 TeV and 5.5 TeV respectively, as obtained 

from a Glauber Monte Carlo, are also given. 
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0.1 Proton-nucleus {pA) collisions 



' 0.11 Glauber formalism 

\C) ' The inelastic cross-section of & p + A reaction, apA, can be derived in the eikonal limit (straight line 

^ . trajectories of colliding nucleons) from the corresponding inelastic nucleon-nucleon AA cross-section, 

^ ! CFm{s) at the center-of-mass energy ^/s, and the geometry of the pA collision simply determined by 

O ' the impact parameter b of the reaction. In the Glauber multiple collision model 1 1 1, such a cross-section 

m ■ 

O 



reads 



]_ _ g-0"Ww('5)TA(6) 
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^ I where Ta(6) is the nuclear thickness function (or nuclear profile function) of the nucleus A at 

Ch ' impact parameter 6: 



>< ■ TA{b) = [ dzpAib,z). 



(2) 



TA{b) gives the number of nucleons in the nucleus A per unit area along a direction z separated 
from the center of the nucleus by an impact parameter b. The nuclear density, pA{b,z), is usually 
parametrized by a Woods-Saxon distribution with nuclear radius Ra = 1-19 • A^/^ — 1.61 • A~^/^ fm 
and surface thickness a = 0.54 fm as given by the experimental data Q and normalized so that 

(fbTA{b) = A. (3) 
0.12 Hard scattering cross-sections 

Though Eq. Q is a general expression for the total inelastic cross-section, it can be applied to an inclu- 
sive p -\- A ^ h -\- X process of production of particle h. When one considers hard scattering processes, 
the corresponding cross-section cr^'^'^ is small and one can expand Eq. Q in orders of a^NTAib) and 
then, to first approximation 
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0.13 "Minimum bias" hard scattering cross-sections 

Integrating Eq. @ over impact parameter, and using one gets the minimum bias {MB) cross-section 
for a given hard process in pA collisions relative to the same cross-section in pp (or A^A^) collisions: 

f ^hard\ a hard /es 

[O-pj^ jMB = A-af^f^ (D) 

From this expression it is easy to see that the corresponding minimum bias multiplicity (invariant 
yield per nuclear reaction: Nj^^^^^ = '^nnpa/^nn.pA^ ^'^^ ^ given hard-process in a.pA collision compared 
to that of a pp collision is 

/ Afhardy a ^NN nrhard ^ hard /f\ 

V^pA /MB — ^ ■ geo ' -'^ NN ~ geo ' "AW ' v"-"/ 
'^pA '^pA 

where a^^" is the geometrical pA cross-section given, in its most general form, by Eq. Q. The 
average nuclear thickness function for minimum bias reactions [making use of Eq. Q] reads: 



/ Ta _ A _ _A 
Thus, for a p-i-Pb {A{Pb) = 208) colhsion at LHC energies ^/s^ = 8.8 TeV with 



J ^ ""pA 



unn ~ 77 mb 13, and (8) 
alll ^ 2162 mb (6l, (9) 

one obtains: {Nppl'^)MB ~ 7.4 • N^ff^, and the average nuclear thickness function amounts to 
{Tpb)MB = 0.096 mb-i = 0.96 fm-^. 



0.2 Nucleus-nucleus {AB) collisions 
0.21 Glauber formalism 

As in the proton-nucleus case, the inclusive inelastic cross-section gab for a collision of nuclei A and B 
is given in the multiple-scattering Glauber approximation by: 



(^ab 



d% 
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-<^NNi^)'^AB('^) 
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where now TAB{b) is the nuclear overlap function of the nuclei A and B separated by impact 
parameter b. TAB{b) can be written as a convolution of the corresponding thickness functions of A and 
B over the element of overlapping area d^s{s = (x, y) is a 2-D vector in the transverse plane, and b is 
the impact parameter between the centers of the nuclei): 



Tab{1>) = j d^sTA{s)TB{\b-s\). 
TAB{b) is normalized so that integrating over all impact parameters one gets: 



(11) 



d^bTAB{b) =AB. 



(12) 



0.22 Hard scattering cross-sections 

As in the pA case, for hard processes of the type A + B ^ h + X, Eq. (flOl . can be approximated by: 



'Tab"^ J d%a'j:,'}[^TAB{b). (13) 



0.23 "Minimum bias" hard scattering cross-sections and yields 

Integrating Eq. (I13t over impact parameter and using (fT2l) . one gets the minimum bias {MB) cross- 
section for a given hard process in AB collisions relative to the corresponding pp cross-section: 



(a^rV= A-B.^-^'' (14) 

'hard 
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Again the corresponding minimum bias multiplicity (invariant yield per nuclear reaction: Nj^"-'^'^ 



'^nnab/^nnab^ ^^'^ ^ given hard-process in a AB collision compared to that of a pp collision is 

I ■\Thard\ /ID ^NN ^jhard ^ ' ^ ^hard n'^\ 

V^AB )mB — a - a ■ — • iV^^ — qeo ' NN ' 

where c7^^° is the geometrical AB cross-section given, in its most general form, by Eq. ( fTUl . The 
average nuclear overlap function for minimum bias reactions [making use of Eq. (fT2t l reads now: 

_ / d'^b Tab A-B AB 

{Tab)mB = f- ,2, - —7^ , p NO = —geo, (16) 

12 



Thus, for a PbH-Pb iA^{Pb) = 43264) collision at LHC energies ^/s^ = 5.5 TeV with 



aNN ~ 72 mb L3J, and (17) 
'^PhPb - 7745 mb (18) 

one gets: {Np^pf)MB ~ 400 • N^'^^'^, and the average nuclear overlap function amounts to 
{TpbPh)MB = 5.58 mh-^ = 55.8 fm'^. 



0.24 Binary collision scaling 

For a given impact parameter b the average hard scattering yield can be obtained by multiplying each 
nucleon in nucleus A against the density it sees along the z direction in nucleus B, then integrated over 
all of nucleus A, i.e. 



(A^aT')(^) = 4^^' / d'^s I pa{s, z') I pB{\b- s\,z") dz"dz' ^ a^Z^ ■ TAB{b) , (19) 

where we have made use of expressions ^ and (fTTT) . In the same way, one can obtain a useful 
expression for the probability of an inelastic NN collision or, equivalently, for the average number of 
binary inelastic collisions, {Ncou), in a nucleus-nucleus reaction with impact parameter b: 

{Ncoii){b) = aMN-TAB{b) (20) 

From this last expression one can see that the nuclear overlap function, TAB{b) = Ncoii{b) / g^n 
can be thought as the luminosity (reaction rate per unit of cross-section) per AB collision at a 



given impact parameter. As an example, the average number of binary collisions in minimum bias Pb+Pb 
reactions at LHC (am = 72 mb = 7.2 fm^) is: 



{NcoIi)mb = 7.2 fm^ • 55.9 fm'^ = 400. (21) 

From (fT9t and d20l . we get so-called "binary (or point-like) scaling" formula for the hard scattering 
yields in heavy-ion reactions: 

(iVj^-^) (6) « {N,oii) (b) • nH^'-'^ (22) 
0.25 Hard scattering yields and cross-sections in a given centrality class 

Equation (IT3l gives the reaction cross-section for a given hard process in AB collisions at a given impact 
parameter 6 as a function of the corresponding reaction cross-section in pp collisions. Very usually, 
however, in nucleus-nucleus collisions we are interested in calculating such a reaction cross-section for 
a given centrality class, (o"^'q'^)ci-C2' where the centrality selection Ci — C2 corresponds to integrating 
Eq. (IT3I between impact parameters hi and 62- It is useful, in this case, to define two parameters |8l|9l: 

• The fraction of the total cross-section for hard processes occurring at impact parameters hi <h < 
62 = 2-Kbdh): 



• The fraction of the geometric cross-section with impact parameter hi < h < b2- 

fgeoibl < 6 < 62) ■ 



27r / \dh(l- e-^NNTABib)] 



, (24) 



[/geo simply corresponds to a O.X (e.g. 0.1) factor for the X%(10%) centrality.] 

Hard scattering production is more enhanced for increasingly central reactions (with larger num- 
ber on Ncoii) as compared to the total reaction cross-section (which includes "soft", - scaling with the 
number of participant nucleons N^art as well as "hard" contributions). The growth with b of the ge- 
ometric cross-section is slower than that of the hard component. For this reason, the behaviour of fj,ard 
and fgeo as a function of b, although similar in shape is not the same (see (8l): fhard ~ 1 for 6 = 2Ra, 
but fgeo ^ 0.75 for b = IRa- 

Similarly to (fT6l . we can obtain now the nuclear overlap function for any given centrality class 
Ci — C2'. 



, _ Ib'dHTAB A-B fhard 

{Tab}c^-C2 = ,2/, = 1?^ ' 1 — 

J^^ d^h <^ab J geo 



The number of hard processes per nuclear collision for reactions with impact parameter hi < b < 
62 is given by 



lAThard\ ^AB {bl < b < 62) . „ (^NN fhard 

^AB [h <b<b2) CJ^s fgeo 



which we could have just obtained directly from ( fT9l and (l25l . From (fT5l and d26l it is also easy 
to see that: 



/ Ajhardx 1 nrhardy fhard 

V^AB )Ci-C2 = V^AB I MB ' —, (27) 

J geo 



Finally, the cross-section for hard processes produced in the centrality class Ci — C2 (correspond- 
ing to a fraction /^eo of the reaction cross-section, (c|g°)ci-C2) is: 

{^'ab'^)ci-C2 = A - B ■ fhard ' cr'mf (28) 

Figure n extracted from 0, plots the (top) fraction of the hard cross-section, fhard{0 < 6 < 62) 
(labeled in the plot as /), as a function of the top fraction of the total geometrical cross-section, /geo(0 < 
6 < 62). for several nucleus-nucleus reactions. 




'^geo 

Fig. 1: Figure 4 of ref. JS^. Fraction of the hard cross-section, f = fhard{0 <b < 62), vs. the fraction of the total geometrical 
cross-section, fgeo{0 <b < bi), for several heavy-ion collisions (from left to right): 197+197, 110+197, 63+193, 27+197, 
and 16+197. 

As a practical application of Eq. d28l and the results of Fig. ^ the hard-scattering cross-sections 
in Pb-i-Pb for the top 0-10% {fhard = 0.41 for fg^o = 0.1 from the practically equivalent Au-i-Au system 
of figure and 0-20% (fhard = 0.664 for fg^o = 0.2) central collisions relate to the pp cross-section, in 
the absence of nuclear effects, respectively as: 

i<^AB\-io% = (208)2. 0.41 • 4^ « 1.7 -10^ -4^ (29) 
i4B\-20% = (208)2-0.664.4^^^ « 2.9-104 -4-^ (30) 



A straightforward way to compute the invariant yield for a given hard process in a given centrality 
class of a nucleus-nucleus collision from the corresponding yield in pp collisions consists in determining, 
via a Glauber MC calculation, the average number of inelastic NN collisions corresponding to that 
centrality class via 

{Ncoll)Ci~C2 = {Tab)ci-C2 ■ ^NN , (31) 



and then use this value in the "binary-scaling" formula 



/ ]\Thard\ 



coll)Ci-C2 ■ ^^NN ; 



or 



(32) 
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{J-AB)Ci~C2- [(^ab)Ci-C2 -^nn 



(33) 



The same two formulae above apply to pA collisions (of course substituting AB by pA and com- 
puting Ncou from Ta instead of from Tab)- 

Finally, to obtain the experimental rates, {J\f^g'''^)ci-C2' actually measured in a given centrality 
bin one needs to take into account the expected integrated luminosity Cint [mb^^] as follows: 



{^flt)c,-C2 



(34) 



0.3 Hard scattering yields and cross-sections for pPh and PhPh collisions at LHC 

As a practical application of the Glauber approach described here, in Tabled the values of {Ncoii) and 
(TpPb) are quoted for different centrality classes obtained from a Monte Carlo calculation flOl for p+Pb 
(y/sNN = 8-8 TeV and a^N = 77 mb) and Pb+Pb (y^SAw = 5.5 TeV for an inelastic pp cross-section of a^^ 
= 72 mb) collisions (Woods-Saxon Pb density parametrization with Ra = 6.78 fm and a = 0.54 fm). 

Table 1: Number of inelastic A^A'' collisions, {Ncoii}, and nuclear thickness (Tppt) or overlap {Tptpb} function per centrality 
class, in p+Pb (y^SAW = 8.8 TeV, ct^W = 77 mb) and Pb+Pb collisions at LHC (^s^w = 5.5 TeV, cr^w = 72 mb) obtained with 
the Glauber Monte Carlo code of ref. 1101 . The errors in (Ncoii), not shown, are of the same order as the current uncertainty in 
the value of the nucleon-nucleon inelastic cross section, crf^f^, at LHC energies (~ 10%). 



Centrahty (Ci - C2) 



p+Pb 







{TpPb) (mb-i) 


{Ncoll) 


{TpbPb) (mb-i) 


0-5% 


15.7 


0.203 


1876.0 


26.0 


0-10% 


15.3 


0.198 


1670.2 


23.2 


10-20% 


13.8 


0.179 


1019.5 


14.2 


20-30% 


12.0 


0.155 


612.4 


8.50 


30-40% 


9.9 


0.128 


351.8 


4.89 


40-50% 


7.8 


0.101 


188.0 


2.61 


50-60% 


5.6 


7.27-10-2 


92.9 


1.29 


60-70% 


3.8 


4.93-10-2 


41.4 


5.75-10-1 


70-80% 


2.6 


3.37-10-2 


16.8 


2.33-10-1 


80-90% 


1.7 


2.20-10-2 


6.7 


9.31-10-2 


90-100% 


1.2 


1.55-10-2 


2.7 


3.75-10-2 


min. bias 


7.4 


9.61-10-2 


400.0 


5.58 



Pb-FPb 



Using d33l . d33l and Table I, we can now easily get the scaling factors of the cross-sections and 
yields from pp to, e.g., central (0-10%), minimum bias, and semi-peripheral (60-80%, from the com- 
bined average 60-70% and 70-80%) p-i-Pb (8.8 TeV) and Pb-i-Pb (5.5 TeV) collisions : 



For p+Pb collisions {(Tm = 2162 mb): 



{N^p"p'-^)o-w% = 15.3 • N^?['' , = 0.198 • 0.1 • 2162 • ^ 450 • 4^ (35) 

(A^,P^'>60-80% = 3.2 • <r , (#)60-80% = 0.042 • 0.2 • 2162 • 4^ « 18 • 4?/^' (36) 

{N';'g)MB = 7.4 • , {a';%';!)MB = 0.096 • 2162 • « 2 • 10^ • a'^t (37) 

For Pb+Pb collisions {all% = 7745 mb): 



(A^S)o-io% = 1670 • <r , i4tpt)o-io% = 23.2 • 0.1 • 7745 • ^ 1-6 • 10^ • *'(38) 
««.)60-80% = 29.1 • iV^^r , i4tPt)m~so% = 0.4 • 0.2 • 7745 • 4^ ^ 6.2 • 10^ • 4^(39) 
{N^tpb)MB = 400 • Nj^^^'' , (4",?i)MiJ = 5.58 • 7745 • 4%' ^ 4.3 • 10^ • (40) 



0.4 Nuclear effects in pA and AB collisions 

Eqs. (|4]i and (fT3b for the hard scattering cross-sections in pj4 and AB collisions have been derived 
within an eikonal framework which only takes into account the geometric aspects of the reactions. Any 
differences of the experimentally measured (Jp'^j^Q ^^^^ respect to these expressions indicate "de facto" 
the existence of "nuclear effects" (such as e.g. "shadowing", "Cronin enhancement", or "parton energy 
loss") not accounted for by the Glauber formalism. Indeed, in the multiple-scattering Glauber model 
each nucleon-nucleon collision is treated incoherently and thus, unaffected by any other scattering tak- 
ing place before (initial-state) or after (final-state effects) it. 



If the Glauber approximation holds, from Q and (fT3l one would expect a. oc A^, and oc A'^ 
growth of the hard processes cross-section with system size respectively. Equivalently, since Nj^^^^g = 
'^nnab/^^naa ^NNAB ~ '^^^^ ~ A^/^, One would expect a growth of the number of hard 
process as oc A^/^ , oc A^/^ for pA^ AA collisions respectively. Experimentally, in minimum bias pA and 
AB collisions, it has been found that the production cross-sections for hard processes actually grow as: 



{cj;T)mb= A'^.a'^f^^ and (^^ )Mfi = • a^' , with a / 1 (41) 

More precisely, in high-pT- processes in pA and heavy-ion collisions at SPS energies one founds 
a > 1 (due to initial-state pj- broadening or "Cronin enhancement"); whereas a < 1 at RHIC ener- 
gies ("high-pT suppression"). Theoretically, one can still make predictions on the hard probe yields in 
pA, AB collisions using the pQCD factorization machinery for the pp cross-section complemented with 
the Glauber formalism while modifying effectively the nuclear PDFs and parton fragmentation functions 
to take into account any initial- and/or final- state nuclear medium effect. 



0.5 Summary of useful formulae 

Finally, let us summarize a few useful formulae derived here to determine the hard-scattering invariant 
yields, cross-sections, or experimental rates, from pp to pA and AB collisions for centrality bin Ci — C2 
(corresponding to a nuclear thickness Ta or nuclear overlap function Tab and to an average number of 
A^A^ inelastic collisions (Ncoii))'- 



/ j2 ]\Thard\ ^2^hard 

-1 -1 - {^A,Ab)Ci-C2 • —j -1— (42) 
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{Ta,Ab)ci-C2 ■ {'^pa.ab)ci-C2 



(]2 hard 
u Upp 
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^int ■ {Ta,Ab)Ci-C2 • {(^pAab)Ci-C2 
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(43) 
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(Q '^pA.ab)Ci-C2 

dprdy 

{d^^^!Zi)c.-C2 
dprdy 



d'^N^I^'''^ 
{Ncou)ci-C2 ■ — -1 — 

dprdy 

,l2 lyrhard 

{Jycoll)Ci-C2 ■ [CrpA.AB)Ci~C2 ■ f^p^f^y 



Cint ■ {Ncou)ci-C2 ■ (f''M°s)Ci-C2 



^2]^hard 

dprdy 



(45) 
(46) 
(47) 
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